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Abstract 

Rationale: During exercise, heart failure patients (HF) show an out-of-proportion ventilation increase, which in patients with 
COPD is blunted. When HF and COPD coexist, the ventilatory response to exercise is unpredictable. 

Objectives: We evaluated a human model of respiratory impairment in 10 COPD-free HF patients and in 10 healthy subjects, 
tested with a progressive workload exercise with different added dead space. We hypothesized that increased serial dead 
space upshifts the VE vs. VC0 2 relationship and that the VE-axis intercept might be an index of dead space ventilation. 

Measurements: All participants performed a cardiopulmonary exercise test with 0, 250 and 500 mL of additional dead 
space. Since DS does not contribute to gas exchange, ventilation relative to dead space is ventilation at VC0 2 = 0, i.e. VE-axis 
intercept. We compared dead space volume, estimated dividing VE-axis intercept by the intercept on respiratory rate axis of 
the respiratory rate vs. VC0 2 relationship with standard method measured DS. 

Main results: In HF, adding dead space increased VE-axis intercept (+0 mL = 4.98 ±1.63 L; +250 mL = 9.69±2.91 L; 
+500 mL= 13.26±3.18 L; p<0.001) and upshifted the VE vs.VC0 2 relationship, with a minor slope rise (+0 mL = 27±4 L; 
+250 = 28±5; +500 = 29±4; p<0.05). In healthy, adding dead space increased VE-axis intercept (+0 mL = 4.9±1.4 L; 
+250 = 9.3±2.4; +500= 13.1 ±3.04; p<0.001) without slope changes. Measured and estimated dead space volumes were 
similar both in HF and healthy subjects. 

Conclusions: VE-axis intercept is related to dead space ventilation and dead space volume can be non-invasively estimated. 
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Introduction 

The behaviour of ventilation during exercise in heart failure 
(HF) and in chronic obstructive pulmonary disease (COPD) 
patients may differ, being characterized in the former by an out-of- 
proportion increase of ventilation (VE), which is greater the 
greater the HF severity [1] and, in the latter, by a normal or 
excessive increase of ventilation in mild or moderate COPD and a 
blunted ventilation increase in severe COPD patients [2-4]. The 
elevated ventilatory response in HF patients seen before lactic 
acidosis ensues and the carbon dioxide (C0 2 ) [5] generated by the 
lactate is trivial relative to the rate of metabolic C0 2 production 
(VC0 2 ) [6,7]. The relationship between VE and VC0 2 is used to 
evaluate ventilatory efficiency [8]; in HF, as well as in pulmonary 
arterial hypertension, an increase of the slope of the VE vs. VC0 2 
relationship is associated with a poor prognosis [9-16]. In COPD, 
ventilatory limitation to exercise is defined either as a reduction of 



ventilatory reserve or as a lowering of inspiratory capacity [17]. In 
case of severe COPD, the rise of ventilation during exercise is 
blunted, and consequently the slope of VE vs. VC0 2 relationship 
is normal or low, being the slope lower the more pronounced the 
emphysema profile [2]. 

HF and COPD often coexist with a reported prevalence of 
COPD in HF patients ranging between 23 and 30% [18] and with 
a relevant impact on mortality and hospitalization rates [19]. In 
patients with COPD and HF, the ventilatory response to exercise 
is poorly predictable. Indeed, HF hyperventilation can be 
counteracted by the incapacity of increasing tidal volume (VT) 
and alveolar ventilation, both being distinctive features of VE 
during exercise in COPD patients [1 7] . As a result, the slope of VE 
vs.VC0 2 relationship might be elevated, normal or even low in 
patients with COPD and HF, regardless of the presence and of the 
severity of ventilatory inefficiency. Up to now, only few studies 
have evaluated the ventilatory behaviour during exercise in 
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patients with coexisting HF and GOPD, being patients with 
comorbidities usually excluded from research trials dedicated to 
HF or COPD [20]. 

In the present study, we evaluated HF patients and healthy 
individuals through a progressive workload exercise with different 
added DS, hoping to mimic at least in part the effects of COPD on 
ventilation behaviour during exercise. We hypothesized that 
increased serial DS upshifts the VE vs. VC0 2 relationship and 
that the VE-axis intercept (VE Y ; nter ) might be an index of DS 
ventilation. Indeed, since DS does not contribute to gas exchange, 
VE relative to DS is VE at VC0 2 = 0, i.e., VE Yint on the VE vs. 
VC0 2 relationship. 

Methods 

Subjects 

Ten HF patients and 10 healthy subjects were enrolled in the 
present study. 

HF patients were regularly followed-up at our HF unit. Study 
inclusion criteria for HF patients were New York Heart 
Association functional classes (NYHA) I to III, echocardiographic 
evidence of reduced left ventricular systolic function (left ventric- 
ular ejection fraction £40%), optimized and individually tailored 
drug treatment, stable clinical conditions for at least 2 months, 
capability/willingness to perform a maximal or near maximal 
cardiopulmonary exercise test (CPET). Patients were excluded if 
they had obstructive and/or restrictive lung disease (forced 
expiratory volume in first second/forced vital capacity ratio 
(FEVi/FVC) <0.70% and/or lung vital capacity (VC) <80% of 
predicted value [21]), clinical history and/or documentation of 
pulmonary embolism, primary valvular heart disease, pulmonary 
artery hypertension, pericardial disease, exercise-induced angina, 
ST changes, severe arrhythmias and significant cerebrovascular, 
renal, hepatic and haematological disease. 

A group of age matched healthy subjects was recruited among 
the hospital staff and from the local community through personal 
contacts. Inclusion criteria were absence of history and/ or clinical 
evidence of any cardiovascular or pulmonary or systemic disease 
contraindicating the test or modifying the functional response to 
exercise, any condition requiring daily medications, and the 
inability to adequately perform the procedures required by the 
protocol. No subjects were involved in physical activities other 
than recreational. 

The investigation was approved by the local ethics committee 
("Ethics committee Centro Cardiologico Fondazione Monzino", 
Institutional Review Board no. S 1 86/ 311) and all participants 
signed a written informed consent before enrolling in the study. 

Study protocol 

At enrolment, demographical and clinical data were collected, 
lung function measurements and echocardiographic evaluation 
were performed to verify that the subjects screened met the study 
inclusion/exclusion criteria, and the informed consent was 
obtained. 

Spirometry (Vmax 29C, SensorMedics, Yorba Linda, CA, US) 
was performed by all participants in accordance with the 
recommended technique [22], and measurements were standard- 
ized as percentages of predicted normal values [23]. 

To become familiar with the procedure, both HF patients and 
healthy subjects had been previously trained to perform an 
exercise test in our laboratory [24]. Thereafter, on different days, 
following a random order, exercise testing was done with 
additional DS equal to 0 mL, 250 mL and 500 mL. 



All participants underwent incremental CPET on an electron- 
ically braked cycle-ergometer (Ergometrics-800, SensorMedics, 
Yorba Linda, CA, US) using a personalized ramp protocol that 
was chosen aiming at a test duration of 10±2 minutes. The 
exercise was preceded by 5 minutes of rest gas exchange 
monitoring and by a 3-minute unloaded warm-up. A 12-lead 
ECG, blood pressure and heart rate were also recorded, and 
arterial oxygen saturation was monitored through a pulse 
oxymeter. The participants wore a nose clip and breathed through 
a mouthpiece connected to a mass flowmeter (Vmax 29C, 
SensorMedics, Yorba Linda, CA, US). Subjects were asked to 
cycle at a pedalling rate of 60-70 rpm, and CPET were self- 
terminated by the subjects when they claimed that maximal effort 
had been achieved. Oxygen consumption (V0 2 ), VC0 2 and VE 
were measured breath by breath with flowmeter and respiratory 
gas sampling lines at the end of the added DS. They were 
averaged every 20 seconds. Anaerobic threshold (AT) was 
calculated with the standard technique [25]. All tests were 
executed and evaluated by 2 expert readers. 

In the absence of psychogenic hyperventilation, below the 
respiratory compensation point [26], the relation between VE and 
VC0 2 is characterized by a linear relationship (VE = aVC0 2 + b), 
with "a" as the slope and "b" as the intercept on the VE axis 
(VE Y int) [8]. Since DS does not contribute to gas exchange, it is 
possible to hypothesize that the ventilation relative to DS is similar 
or related to the VE at VC0 2 = 0, which is the Y intercept of VE 
vs. VC0 2 relationship. To calculate DS volume (VD) from VE Y i n t 
(VD Y ; nt ), we need to identify the corresponding respiratory rate 
(RR). This was obtained as the intercept of the RR vs. VC0 2 
relationship on the RR axis (RR Y ; nt ). Specifically, the RR vs. 
VC0 2 relationship was calculated through its linear portion that 
starts from the beginning of exercise and ends when RR increases 
more steeply, which corresponds to the tidal volume inflection/ 
plateau [27,28]. An example on how we calculate VE Y i n t and 
RR Y i nt is reported in figure 1 . 

We compared estimated VD values (VD Yint ) with resting and 
exercise values of VD, measured with standard method [8] 
(VD meas ), in the 3 experimental conditions, with 0 mL, 250 mL 
and 500 mL of added DS. The volume of mouthpiece and 
flowmeter (50 mL) was subtracted from VD. The standard 
calculation of VD [8] (VD meas ) is obtained by the following 
equation: 

VD = VT*[1 - -(863 * VC02)/(VE* PaC02)] 

with 863 as a constant and P a C0 2 as pressure for arterial C0 2 . 

In healthy individuals [29], but not in HF patients [30], P a C0 2 
can be reliably estimated from end-tidal expiratory pressure for 
C0 2 (P ET C0 2 ). Therefore, we measured P a C0 2 from arterial gas 
sampling in HF patients, and we estimated P a C0 2 from P EX C0 2 
in healthy subjects. Thus, only in HF patients, a small catheter was 
introduced into a radial artery, blood samples were obtained at 
rest and every 2 minutes during exercise, and P a C0 2 was 
determined with a pH/blood gas analyzer (GEM 4000, Instru- 
mentation Laboratory, Bedford, MA, US). 

We calculated possible VD changes during exercise, and we 
evaluated whether an added DS modifies the slope of the VE vs. 
VC0 2 relationship and/or it simply upshifts it. 

Statistical analysis 

Data are mean ± standard deviation (SD). Cardiopulmonary 
measurements were collected breath by breath and reported as 
average over 20 s. Comparisons between the two groups were 
done through unpaired t-test. Both in HF and in healthy subjects, 
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Figure 1 . VE vs. VC0 2 relationship in a patient. The relationship is 
linear up to the respiratory compensation point (end of the isocapnic 
buffering period) (Upper panel). RR vs. VC02 relationship. The 
relationship is calculated as for VE vs. VC0 2 (Lower panel). 
doi:1 0.1 371 /journal.pone.0087395.g001 

analysis of variance for repeated measures with Bonferroni post 
hoc test was performed to analyze the effect of the adding of 
different DS and to evaluate the changes of VD meas during 
exercise in the 3 experimental conditions. Bland and Altaian 
relationship was calculated to compare VD Y ; m values and VD meas 
values in HF patients and in healthy individuals. 

Statistical significance was set at p<0.05. All statistics were 
performed with IBM SPSS statistics 20.0 for windows. 

Results 

We enrolled 1 0 HF patients (9 males; mean age 6 1 ± 1 3 years) 
and 10 age-matched healthy subjects (8 males; mean age 59± 10 
years). The main anthropometric data were not significandy 
different between the two groups. Patients with HF and healthy 
subjects were free from obstructive defects; although within the 
predicted normal limits, lung volumes tended to be smaller in HF 
patients than in normal subjects (table 1). 

HF patients 

Mean left ventricle ejection fraction was 33±5%. The cause of 
HF was ischemic dilated cardiomyopathy in 4 cases and primary 
dilated cardiomyopathy in 6 cases. Three patients had an 
implantable cardioverter defibrillator; 9 were in sinus rhythm 
and 1 was in permanent atrial fibrillation. Four patients were in 



NYHA class I, 5 in NYHA class II and 1 in NYHA class III. All 
HF patients were on (3-blockers, 9 with angiotensin-converting 
enzyme inhibitors, 4 with aldosterone receptor antagonists, 5 with 
diuretics and 3 with amiodarone. 

All HF patients performed CPET without added DS and with 
250 mL and 500 mL of additional DS without complications. In 
the HF group, peak VO2 was slightly reduced compared to 
healthy subjects. With the exception of reduced peak workload 
and of an increased VT, the adding of different DS did not 
significandy impact on CPET data at peak of exercise and on V0 2 
at AT (table 2). In table 3 VE, RR, VT, VD/VT, VC0 2 , P ET C0 2 
and P a C0 2 during exercise are reported with 0, 250 and 500 mL 
of added DS. 

Values of VEy int , RR Yint , VD Yint , VD meas and the slope of VE 
vs VC0 2 relationship in HF patients with 0 mL, 250 mL and 
500 mL of additional DS are reported in table 4. 

With the adding of DS, the VE Y i„ t increased significantly, 
whereas RRyint showed a limited increase. Adding DS upshifted 
the VE vs. VC0 2 relationship with a minor slope increase 
(figure 2). 

The calculated VD Y ; nt rose as added DS increased; mean 
VD Y mt increase with 250 and 500 mL of added space was 
226±127mL and 446±123mL. VD meas increased during 
exercise in the 3 conditions albeit only as a trend when DS was 
not added (table 5). 

Figure 3 reports the Bland and Altaian plot of VD Y i nt vs. 
VD meas at rest for HF patients in the 3 exercise conditions. As an 
average, a good agreement was observed when VD was calculated 
either by VE Y i nt , or VD meas , with or without additional DS. 

Healthy subjects 

Healthy subjects performed all CPET without complications. 
Peak exercise data and V0 2 at AT were not significantly affected 
by the adding of DS (table 2). 

When DS was added, the value of the slope of VE vs. VC0 2 
relationship and RRyint did not change, whereas only the VE Yint 
increased significandy (table 4) with an upshift of the relationship 
(figure 4). Similarly to HF patients, VD Yint increased with added 
DS in the three experimental conditions, specifically by 
300±150mL and by 570±160mL with 250 and 500 mL, 
respectively. 

During exercise, VD meas remained constant without additional 
DS, whereas it significantly decreased during exercise with added 
DS, but this finding is likely due to the underestimation of P a C0 2 
by P ET C0 2 with added DS (table 5). 

Figure 5 reports the Bland and Altaian plot of VD Y i„ t vs. 
VD meas at rest for healthy subjects and showed a good correlation 
between the two methods both with and without additional DS. 

Discussion 

In the present study, we evaluated a human model of increased 
dead space in HF patients and in healthy subjects, applying a 
progressive workload exercise with different added DS. We 
documented that a rise in serial DS, mimicking a rise in 
anatomical DS, was parallel to the VE Yint increase both in healthy 
individuals and in HF patients. Therefore, VE Y ; nt is related to DS 
ventilation. Moreover, we showed that the value of DS can be 
non-invasively estimated as the ratio of VE Y ; nt /RR Y ; nt . 

Few study limitations should be discussed at first. Firstly, our 
research was undertaken to analyze the role on ventilation 
behaviour during exercise of a respiratory comorbidity, COPD, 
in HF patients. We built a COPD model by adding an external 
dead space. We recognize that our model is only a partial COPD 
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Table 1. Main anthropometric characteristics, demographical and pulmonary function data of heart failure patients and healthy 
subjects enrolled in the study. 







HEART FAILURE PATIENTS 


HEALTHY SUBJECTS 


p value 


Number 


10 


10 


NS 


Male/female 


9/1 


8/2 


NS 


Age (yr) 


61 ±12 


59±7 


NS 


Height (cm) 


172±9 


173±6 


NS 


Weight (Kg) 


85±15 


77±11 


NS 


BMI (Kg/m 2 ) 


28.6±3.8 


25.4±3.2 


NS 


VC (L) 


3.58±0.75 


4.72 ±1.03 


<0.01 


VC (% predicted) 


91 ±14 


112±13 


<0.01 


FVC (L) 


3.47 ±0.67 


4.63±1.10 


<0.01 


FVC (% predicted) 


90±12 


112±14 


<0.01 


FEV, (L) 


2.56±0.58 


3.57 ±0.84 


<0.001 


FEVi (% predicted) 


79±14 


107±17 


<0.001 


FEV, /FVC 


73 ±4 


76±5 


NS 



Data are presented as number or mean ± SD. BMI = body mass index; NS = not significant; FEV, = forced expiratory volume in 1 s; FVC = forced vital capacity; VC = 
vital capacity. 

doi:1 0.1 371 /journal.pone.0087395.t001 



model because we have not considered any of the systemic 
consequences of COPD and we have limited our attention to DS 
changes. Our model was over-simplistic also as regards lung 
mechanics because an artificial dead space increase does not 



generate air trapping which is one of the most characteristic 
features of COPD during exercise. Secondly, our model was short 
lasting, so that chronic ventilatory and chemoreceptor adaptations 
to increased DS were not evaluated as were not evaluated 



Table 2. Cardiopulmonary exercise testing data in heart failure patients (upper panel) and healthy subjects (lower panel) with 
0 mL, 250 mL and 500 mL of additional dead space. 





HEART FAILURE PATIENTS 


ADDED DEAD SPACE 






ANOVA p value 




+0 mL 


+250 mL 


+500 mL 




Peak workload (W) 


109±41* 


103±47 


96±41 


0.006 


Peak V0 2 (mL/min/Kg) 


19.9±5.8 


19.3 ±5.6 


19.6±5 


NS 


V0 2 at AT (mL/min/Kg) 


13±3 


14.1 ±4 


12.7±5.8 


NS 


Peak 0 2 pulse (mL/beat) 


15.8±5.7 


15.4±5.2 


15.7±4.8 


NS 


Peak HR (beat/min) 


1 1 1 ±26 


110±28 


104±20 


NS 


Peak VT (L) 


1.9±0.49 


1.93±0.49 $ 


2.09±0.59 


0.047 


Peak VE (L/min) 


55.6±14 


59.8±14 


58.8±11 


NS 


Peak RR (bpm) 


30±4 


31 ±5 


30±5 


NS 


Peak Pa0 2 (mmHg) 


107±12 


104±16 


100±20 


NS 


Peak Sa0 2 (L/min) 


98.4±1.2 


97.5 ±1.9 


97.7±1.7 


NS 


HEALTHY SUBJECTS 


Peak workload (W) 


200±51 


195±51 


189±45 


NS 


Peak V0 2 (mL/min/Kg) 


36.1 ±8.4 


35.6±7.2 


35.8±7.5 


NS 


V0 2 at AT (mL/min/Kg) 


21.7±5.7 


23.6±3.7 


25.3±6.6 


NS 


Peak 0 2 pulse (mL/beat) 


17.5±4.2 


17±2.9 


18.4±3.4 


NS 


Peak HR (beat/min) 


156±18 


157±18 


156±18 


NS 


Peak VT (L) 


2.71 ±0.6 


2.57 ±0.9 


2.95±0.5 


NS 


Peak VE (L/min) 


88.6±21.9 


87.2±16.2 


88.6±17.1 


NS 


Peak RR (bpm) 


32±4 


32±6 


30±5 


NS 



Data are presented as means ± SD; AT = anaerobic threshold; bpm - breaths per minute; HR = heart rate; NS= not significant; P a 0 2 = arterial oxygen pressure; 
RR = respiratory rate; S a 0 2 = arterial oxygen saturation; RR = respiratory rate; VC* 2 = oxygen consumption; VE = ventilation; VT = tidal volume; W = watt. 
*p<0.05 versus +500 mL; * p<0.01 versus +500 mL. 
doi:1 0.1 371 /journal.pone.0087395.t002 
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Table 3. Ventilatory parameters in heart failure patients with 0, 250 and 500 mL of additional dead space. 



HF PATIENTS 


+0 mL 




+250 mL 


+500 mL 


ANOVA p value 


Rest 


VE (Umin) 


11.8 




1.7*" 


16.2 




3.5 


20.0 




4.2 


<0.001 


RR (bpm) 


14.2 




2.0 


16.4 




4.1 


16.8 




3.1 


NS 


VT(L) 


0.8 




0.2* 


1.0 




0.2 £ 


1.2 




0.1 


<0.001 


VD/VT 


0.47 




0.1 5 $& 


0.61 




0.10 


0.67 




0.11 


< 0.001 


VC0 2 (LVmin) 


0.25 




0.06 


0.29 




0.13 


0.29 




0.14 


NS 


PetC0 2 (mmHg) 


33.4 




1.6 


33.0 




2.5 


33.1 




4.2 


NS 


P a C0 2 (mmHg) 


35.8 




22 $n 


38.6 




1.9 


39.9 




2.02 


<0.001 


4 min exercise 


VE (LVmin) 


21.6 




3.8u # 


30.2 




5.0 


34.8 




4.3 


<0.001 


RR (bpm) 


18.7 




2.7 


20.4 




4.3 


20.7 




4.1 


NS 


VT(L) 


1.2 




0.2 s 


1.5 




0.3 


1.7 




0.3 


<0.001 


VD/VT 


0.33 




0.09 $ " 


0.45 




0.06 


0.54 




0.10 


< 0.001 


VC0 2 (LVmin) 


0.64 




0.15 


0.74 




0.17 


0.73 




0.21 


NS 


PetC0 2 (mmHg) 


37.2 




2.9 


35.7 




3.6 


37.4 




4.2 


NS 


P a C0 2 (mmHg) 


38.4 




2.8 


38.8 




3.4 


41.2 




3.9 


NS 


8 min exercise 


VE (LVmin) 


39.9 




5.9" 


44.5 




4.8 £ 


52.4 




8.4 


<0.001 


RR (bpm) 


25.1 




3.2 


25.3 




5.2 


26.8 




4.6 


NS 


VT(L) 


1.6 




0.3 


1.8 




0.4 


2.0 




0.5 


NS 


VD/VT 


0.28 




0.0& 1 * 


0.41 




0.07 


0.46 




0.09 


< 0.001 


VC0 2 (Umin) 


1.28 




0.35 


1.27 




0.29 


1.34 




0.35 


NS 


PetC0 2 (mmHg) 


37.2 




4.3 


36.8 




4.6 


38.5 




4.2 


NS 


P a C0 2 (mmHg) 


38.0 




3.7 


39.4 




4.2 


41.4 




4.6 


NS 


peak exercise 


VE (Umin) 


55.7 




14.0 


59.9 




14.6 


58.9 




11.3 


NS 


RR (bpm) 


30.3 




4.7 


31.4 




4.0 


29.8 




5.0 


NS 


VT(L) 


1.9 




0.5 


1.9 




0.5 


2.1 




0.6 


NS 


VD/VT 


0.26 




0.1 1»" 


0.39 




0.10 


0.45 




0.11 


< 0.001 


VC0 2 (Umin) 


1.81 




0.67 


1.72 




0.68 


1.58 




0.55 


NS 


PetC0 2 (mmHg) 


35.4 




4.5 


35.64 




4.8 


39.0 




4.9 


NS 


P a C0 2 (mmHg) 


35.8 




3.8 


38.0 




4.2 


41.3 




5.5 


0.049 



Data are presented as means ± SD; VE = ventilation; RR = respiratory rate; VT = tidal volume; VD = dead space volume; VC0 2 = 

arterial carbon dioxide pressure; P ET C0 2 = End-tidal carbon dioxide pressure; bpm- breaths per minute; 

$ : p<0.05 vs. 250 mL; ": p<0.001 vs. 500 mL; *: p<0.001 vs. 250 mL; & : p<0.01 vs.500 mL; # <0.01 vs. 250 mL. 

doi:1 0.1 371 /journal.pone.0087395.t003 



carbon dioxide production; P a C0 2 = 



primitive chemoreceptor abnormalities as drivers of the alveolar 
hypoventilation observed in COPD patients. Thirdly, with the Y- 
intercept we analyze an index of overall DS. However, in the 
present setting, we were able to change DS only by adding an 
external (anatomical equivalent) DS, so that we do not know if 
changes in physiological DS similarly influence the VEyint- 
Fourthly, VE changes during exercise are due to VC0 2 , VD/ 
VT and P a C0 2 changes, and all may influence the VE vs. VC0 2 
relationship. In the present study, we added external DS, which at 
each step of exercise, was associated to an increase of VD/VT and 
P a C0 2 (the latter in 2 steps only as a trend) resembling what 
happens during exercise in COPD patients (table 3). Therefore 
both P a C0 2 and VD/VT changes have likely a role in the VE vs. 
VC0 2 relationship changes we observed after adding DS. It is 
recognized that P a C0 2 measurements were done only in HF 
patients and not in healthy subjects, but a different behaviour in 



healthy subjects is unlikely. Fifthly, the condition of VE at C0 2 
production equal 0, as such at the VE Y mt of the VE vs. VG0 2 
relationship, is a mathematical extrapolation with no physiological 
meaning. Moreover, absolute DS changes during exercise, so that 
also the VE Y mt value is likely close but different from the rest 
value. Indeed, we showed that VD tended to increase in HF 
patients and to reduce in healthy subjects during exercise without 
added DS. However, we suggest using VE Y in t as a tool to evaluate 
the presence of an increased DS, regardless of its physiological 
meaning with respect to rest and exercise. The adding of DS 
significantly reduced the external work produced in HF patients, 
while a not significant reduction was observed in normal subjects. 
Peak V0 2 remained unchanged in both groups after adding DS; 
this finding suggests that added DS was associated to an increased 
work of breathing which, as a percentage of total work, seems to be 
greater in HF patients than in normal subjects. 
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Table 4. Values of the slope of VE vs VC0 2 relationship, VE Yin t, RRyint ar| d volume of dead space in heart failure patients (upper 
panel) and healthy subjects (lower panel) with 0 mL, 250 mL and 500 mL of additional dead space. 





HEART FAILURE PATIENTS 


ADDED DEAD SPACE 






ANOVA p value 




+0 mL 


+250 mL 


+500 mL 




VE/VC0 2 slope 


27±4 


28±5 


29±4 


0.037 


VE Yin , (L/min) 


4.98±1.63 t5 


9.69±2.91* 


13.26±3.18 


0.000 


RR Yint (bpm) 


13±4 &$ 


15+3 


16+3 


0.032 


VD Yin , (L) 


0.39±0.07 s>s 


0.61 ±0.1 2 s 


0.83 ±0.11 


0.000 


VD meas (L) 


0.38±0.08 OS 


0.61 ±0.1 2 s 


0.80±0.09 


0.000 


HEALTHY SUBJECTS 


VE/VC0 2 slope 


23±3 


24±4 


24+4 


NS 


VE Yin , (L/min) 


4.9±1.4 ts 


9.3 ± 2.4 s 


13.1 ±3.04 


0.000 


RR Yint (bpm) 


14±4 


14±4 


14±3 


NS 


VD Yin , (L) 


0.37 ±0.11 05 


0.68±0.15 s 


0.95 ±0.1 4 


0.000 


VD meas (L) 


0.37±0.06 OS 


0.68±0.11* 


0.94 ±0.1 


0.000 


Data are presented as means ± SD; RR Yint = respiratory rate calculated as Y intercept of RR vs VC0 2 relationship; VC0 2 = carbon dioxide production; VD Yint = dead 
space volume calculated as VE Yint /RR Y jnt; VD meas = dead space volume measured by P a C0 2 in heart failure patients and estimated by PfrC0 2 in healthy subjects; VE = 
ventilation; VE Y j nt = ventilation at VCO 2 = 0, calculated as Y intercept of VE vs VC0 2 relationship. 
t p<0.001 versus +250 mL; 
5 p<0.001 versus +500 mL; 



*p<0.01 versus +500 mL; 
& p<0.05 versus +250 mL; 
$ p<0.05 versus +500 mL; 
°p<0.01 versus +250 mL. 
doi:l 0.1 371 /joumal.pone.0087395.t004 



VE Yint (L/min) i; 
13.26±3.1S 

10 

9.69±2.91* 
4.98±1.63 



29±4= SLOPE 



28*5" SLOPE 



27±4 = SLOPE 




::■ :■ 



400 



600 



S00 



1000 



1200 
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Figure 2. VE vs. VC0 2 relationship in heart failure patients with 0 mL (black line), 250 mL (grey line) and 500 mL (dotted line) of 
additional dead space (DS). The adding of DS uplifts the VE vs. VC0 2 relationship with a minor slope increase, f p<0.001 versus +250 mL; § 
p<0.001 versus +500 mL; * p<0.01 versus +500 mL; * p<0.05 versus other all. 
doi:1 0.1 371 /journal.pone.0087395.g002 
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Table 5. Values of volume of dead space at rest and during exercise in heart failure patients and healthy subjects with no 
additional dead space and with 250 mL and 500 mL of additional dead space. 









+0 n 


iL 










+250 


mL 






+500 mL 










HF 






H 






HF 




H 




HF 


H 




VD 

rest (L) 




0.38 




0.08 


0.37 




0.06 


0.61° 


± 0.12 


0.68 ¥ ^ 


± 0.1 1 


0.80 12 ± 0.09 


0.94 pWS 


£ °± 0.10 


VD meas 2 


(L) 


0.38 




0.11 


0.36 




0.04 


0.63 " 


± 0.07 


0.57 ± 


0.13 


0.87 ± 0.08 


0.70 ± 


0.17 


VD meas 4 


(L) 


0.39 




0.08 


0.34 




0.05 


0.68 " 


± 0.11 


0.56 ± 


0.09 


0.91 ± 0.09 


0.67 ± 


0.16 


VD meas 6 


(L) 


0.43 




0.19 


0.36 




0.08 


0.71 


± 0.13 


0.51 ± 


0.09 


0.92 ± 0.15 


0.62 ± 


0.15 


VD meas 8 


(L) 


0.43 




0.09 


0.32 




0.08 


0.73 " 


± 0.11 


0.48 ± 


0.12 


0.90 ± 0.14 


0.57 ± 


0.12 


VD meas 
peak (L) 




0.45 




0.18 


0.31 




0.11 


0.71 " 


± 0.13 


0.44 ± 


0.08 


0.90 ± 0.13 


0.55 ± 


0.12 


p value 




NS 






NS 






0.001 




0.001 




0.05 


0.001 





Data are presented as means ± SD; DS = dead space; H = healthy subjects; HF = heart failure patients; NS = not significant; VD Y i„ t = dead space volume calculated as 

VE Y jnt/RRYint; VDmeas = dead space volume measured by P a C0 2 in heart failure patients and estimated by P ET C0 2 in healthy subjects. 

°p<0.001 versus VD m(!as 6'; Q p<0.05 versus VD m< , as 6'; ¥ p<0.05 versus VD meas 6'; ,l p< 0.001 versus VD meas 8'; • p<0.01 

versus VD meas peak; p p<0.001 versus VD meas 2'; °° p<0.001 versus VD meas 4'; € p<0.001 versus VD meas 6'; 

£ p<0.001 versus VD meas 8'; " p<0.001 versus VD meas peak. 

doi:1 0.1 371/journal.pone.0087395.t005 



We measured DS during exercise using a standard formula [8] 
in HF patients. To avoid systemic artery catheterization, we 
estimated P a C0 2 from PetC0 2 in healthy subjects, which is a 
accepted method in the absence of lung disease [29]. It is 
recognized, however, that albeit largely used in the clinical setting, 
extrapolation of P a C0 2 from P EX C0 2 even in normal individual is 
approximate and likely to cause some of the variability observed 
(figure 5). Moreover, the values obtained in normal subjects with 
added DS showed a progressive and unrealistic DS reduction. This 
is due to a P a C0 2 underestimation by PetC0 2 when adding DS, 
confirming the need to directly measure P a C0 2 during exercise for 
DS evaluation [30]. The low PetC0 2 compared to P a C0 2 
observed during exercise with added dead space (Table 3) is likely 



due to the rapid rise of PC0 2 during exhalation, which does not 
reach a plateau. 

Adding DS increased the slope of VE vs. VC0 2 relationship in 
HF patients but not in control subjects. This is different from what 
happens in patients with severe COPD who show a high VE/ 
VC0 2 ratio at the beginning of exercise but a blunted VE increase 
during exercise, so that the slope of VE vs. VG0 2 relationship is 
normal or low [2] . In our model, the DS increase was too modest 
to generate a ventilatory limitation to exercise, being the 
ventilatory reserve at peak exercise always preserved. Accordingly, 
in HF patients, but not in healthy subjects, we observed a minor 
exercise performance reduction with the adding of DS. 

The VE vs. VC0 2 relationship is frequently used as a 
prognostic tool in HF patients [9-13]. Some laboratories prefer 
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~ 0.3 
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Figure 3. Bland and Altman plot of estimated dead space (DS) volume calculated as VE Yinl /RRYint |VD vin ,| and measured DS volume 
(VD mea s) at rest, calculated as (1-863/P a C0 2 (VE/VC0 2 ) x VT) for heart failure patients with 0 mL (diamonds), 250 mL (circles) and 
500 mL (crosses) of additional DS. The grey line identifies the mean difference of VD meas - VD Yi nt; the black lines identify the mean difference of 
VD meas - and VD Y i n t± 1 .96*standard deviation. P a C0 2 = arterial carbon dioxide pressure; VE = ventilation; VT = tidal volume. 
doi:1 0.1 371 /journal. pone.0087395.g003 
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24±4 SLOPE 



24±4 SLOPE 



23±3 SLOPE 




VE Ylllt (L'min) 

15 



VCO : (mL) 

Figure 4. VE vs. VC0 2 relationship in healthy subjects with 0 mL (black line), 250 mL (grey line) and 500 mL (dotted line) of 
additional dead space (DS). The adding of DS upshifts the VE vs VC0 2 relationship without significant slope changes, f p<0.001 versus +250 mL; 
§ p<0.001 versus +500 mL. 
doi:1 0.1 371 /journal.pone.0087395.g004 



to analyze the ratio of the relationship [31], others the slope [32]. 
However, the ratio varies during exercise, so that which exercise 
VE/VC0 2 ratio value should be considered is still a matter of 
debate [31]. Moreover, while the behaviour of VE/VCO z ratio 
during exercise is well described in normal and HF individuals 
[31], its behaviour in COPD or in patients with HF and COPD is 
less characteristic and not used as a diagnostics/prognostic tool. 
To avoid the above-mentioned uncertainties, many authors prefer 
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Figure 5. Bland and Altman plot of estimated dead space (DS) 
volume calculated as VE Yint /RR Y rnt (VD Yin ,) and measured DS 
volume (VD meas ) at rest, calculated as (1 -863/P a C0 2 (VE/ 
VC0 2 )*VT) with P a C0 2 for healthy subjects with 0 mL (dia- 
monds), 250 mL (circles) and 500 mL (crosses) of additional 
DS. The grey line identifies the mean difference of VD meas - VD Y i n t ; the 
black lines identify the mean difference of VD meas - and VD Yint ±1.96*- 
standard deviation. P a C0 2 was estimated from P ET C0 2 . P a C0 2 = carbon 
dioxide pressure; P ET C0 2 = tele-expiratory carbon dioxide pressure; 
VE = ventilation; VT = tidal volume. 
doi:1 0.1 371 /journal.pone.0087395.g005 



to study the VE vs. VC0 2 relationship throughout the exercise 
[33] or up to the respiratory compensation point [8], To do so, the 
slope of the VE vs. VC0 2 relationship is calculated, but no 
attention is dedicated to the intercept of this relationship on the 
VE axis. However, the increase of the slope of VE vs. VC0 2 
relationship may be blunted when COPD is associated to HF [2]. 
Notably, the presence of COPD in HF may be difficult to be 
defined because some lung impairment is typical of HF and 
particularly in more advanced cases regardless of COPD [5]. In 
the present study, we showed that a DS increase is parallel to the 
VE Y int increase, so that its value should be taken into account 
when analyzing the VE vs. VC0 2 relationship. Indeed, VEy mt 
differences were observed even by adding a relatively small DS 
(250 mL), which corresponded to 1/10 of peak VT in healthy 
subjects. It is recognized, however, that whilst the means of 
estimated and measured VD are similar, the individual values 
differ up to 60% in case of no added DS and up to ~20% when 
500 mL DS were added. This suggests caution when analyzing 
specific individual data, particularly in the presence of no or 
modest lung disease. 

In the present study, we added 250 mL and 500 mL of DS 
during exercise. To confirm that VEymt increase was related to DS 
increase, we calculated VD Ymt . To do so, we need to divide VE by 
RR, but the value of RR to be chosen is an open question. We 
used the intercept of the RR vs. VC0 2 relationship on the RR axis 
because this is the RR value corresponding to VEy int . Interest- 
ingly, the changes of VD Ymt values with added DS were very 
similar to the amount of added DS. 

In conclusion, we provide the rational basis for the assessment of 
VEyint during exercise as a tool to evaluate DS. Further studies are 
needed to confirm and to analyze the clinical meaning of the 
present observation. 
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At a Glance Commentary 

The ventilation (VE) vs. VCO z relationship during exercise is 
commonly used to assess ventilatory efficiency and prognosis in 
heart failure patients. The slope of the VE vs. VC0 2 relationship 
increases as heart failure severity increases, whereas in respiratory 
patients the VE vs. VC0 2 slope during exercise is reduced the 
greater the ventilatory limitation. However, respiratory disease 
often coexists in heart failure patients so that the mean of the slope 
of the VE vs. VC0 2 relationship in these cases is unclear. 

We reasoned that the VE vs. VC0 2 behavior during exercise is 
a linear relationship, at least up to the respiratory compensation 
point, characterized by a slope and a Y intercept value. The latter 
has been ignored, but it represent the ventilation at VOC 2 = 0 and 
therefore it is somehow related to dead space ventilation. 
Accordingly, we built a human model of increased anatomical 
dead space, resembling what happens in chronic obstructive 
pulmonary disease, by adding external dead space during exercise 
in healthy subjects and HF patients. We demonstrated that adding 

References 

1. Metra M, Dei Cas L (1996) Role of cxereise ventilation in the limitation of 
functional eapaeity in patients with congestive heart failure. Basie Res Cardiol 91 
Suppl 1: 31-36. 

2. Paoletti P, Dc Filippis F, Fraioli F, Cinquanta A, Valli G, et al. (2011) 
Cardiopulmonary exercise testing (CPET) in pulmonary emphysema. Respir 
Physiol Neurobiol 179: 167-173. 

3. Pierce AK, Luterman D, LoudermilkJ, Blomqvist G, Johnson RL Jr (1968) 
Exercise ventilatory patterns in normal subjects and patientws with airway 
obstruction. J Appl Physiol 25: 249-254. 

4. Teopompi E, Tzani P, Aiello M, Ramponi S, Visca D, et al. (2013) Ventilatory 
Response to Carbon Dioxide Output in Patients with Chronic Heart F ailure and 
in Patients with Chronic Obstructive Pulmonary Disease with Comparable 
Exercise Capacity. Respir Care. 

5. Wasscrman K, Zhang YY, Gitt A, Belardinelli R, Koike A, et al. (1997) Lung 
function and exercise gas exchange in chronic heart failure. Circulation 96: 
2221-2227. 

6. Chua TP, Harrington D, Ponikowski P, Webb-Pcploc K, Poole- Wilson PA, et al. 
(1997) Effects of dihydrocodeine on chemosensitivity and exercise tolerance in 
patients with chronic heart failure. J Am Coll Cardiol 29: 147—152. 

7. Picpoli M, Clark AL, Volterrani M, Adamopoulos S, Sleight P, et al. (1996) 
Contribution of muscle afferents to the hemodynamic, autonomic, and 
ventilatory responses to exercise in patients with chronic heart failure: effects 
of physical training. Circulation 93: 940-952. 

8. Wasscrman K, Hansen JE, Sue DY, Stringer WW, Sictscma KE, et al. 
(Principles of exercise testing and interpretation. 2012. 5th ed. Lippincott 
williams & wilkins.). 

9. Arena R, Myers J, Abella J, Peberdy MA, Bensimhon D, et al. (2007) 
Development of a ventilatory classification system in patients with heart failure. 
Circulation 115: 2410-2417. 

10. Corra U, Mezzani A, Bosimini E, Scapellato F, Imparato A, et al. (2002) 
Ventilatory response to exercise improves risk stratification in patients with 
chronic heart failure and intermediate functional capacity. Am Heart J 143: 
418-426. 

11. Francis DP, Shamim W, Davics LC, Piepoli MF, Ponikowski P, et al. (2000) 
Cardiopulmonary exercise testing for prognosis in chronic heart failure: 
continuous and independent prognostic value from VE/VCO(2)slope and peak 
VO(2). Eur Heart J 21: 154-161. 

12. Gitt AK, Wasserman K, Kilkowski C, Klecmann T, Kilkowski A, et al. (2002) 
Exercise anaerobic threshold and ventilatory efficiency identify heart failure 
patients for high risk of early death. Circulation 106: 3079-3084. 

1.3. Richer FX, Victzkc G, Wcrneckc KD, Bauer U, Opitz C, ct al. (2000) 
Impairment of ventilatory efficiency in heart failure: prognostic impact. 
Circulation 101: 2803-2809. 

14. Paolillo S, l r arina S, Bussotti M, Iorio A, PerroncFilardi P, et al. (2012) Exercise 
testing in the clinical management of patients affected by pulmonary arterial 
hypertension. Eur J Prev Cardiol 19: 960-971. 

15. Deboeck G, Scoditti C, Huez S, VachicryJL, Lamottc M, ct al. (2012) Exercise 
testing to predict outcome in idiopathic versus associated pulmonary arterial 
hypertension. Eur Respir J 40: 1410-1419. 

16. Schwaiblmair M, Faul C, von Scheidt W, Berghaus TM (2012) Ventilatory 
efficiency testing as prognostic value in patients with pulmonary hypertension. 
BMC PulmMcd 12: 2.3. 



dead space increases the Y intercept of the VE vs. VC0 2 
relationship. The Y intercept of VE vs. VC0 2 relationship is 
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